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Abstract 

The proper way of averaging is an important question with regards to Tsal- 
lis' Thermostatistics. Three different procedures have been thus far employed 
in the pertinent literature. The third one, i.e., the Tsallis-Mendes-Plastino 
(TMP) Q normalization procedure, exhibits clear advantages with respect to 
earlier ones. In this work, we advance a distinct (from the TMP-one) way of 
handling the Lagrange multipliers involved in the extremization process that 
leads to Tsallis' statistical operator. It is seen that the new approach consid- 
erably simplifies the pertinent analysis without losing the beautiful properties 
of the Tsallis-Mendes-Plastino formalism. 
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I. INTRODUCTION 



Tsallis' thermostatistics [|l|-|6| is by now recognized as a new paradigm for statistical 
mechanical considerations. One of its crucial ingredients, Tsallis' normalized probability 
distribution PJ, is obtained by following the well known MaxEnt route |7J. One maximizes 
Tsallis' generalized entropy 0.|[ 

k ~ q-l ' 1 ' 

(k = k(q) tends to the Boltzmann constant k^ in the limit q — > 1 f2|) subject to the 
constraints (generalized expectation values) 

w 

£ft = 1 ( 2 ) 

i=l 

=«Qi» g , (3) 

where is the probability assigned to the microscopic configuration i (i — 1, . . . , w) and one 
sums over all possible configurations w. Oj (j = 1, . . . , n) denote the n relevant observables 
(the observation level 0), whose generalized expectation values ((Oj)) are (assumedly) a 
priori known. 

The Lagrange multipliers recipe entails maximizing [|TJ 

C / w \ n ( y^to rfln^> \ 



yielding 



where 



ji/(i-g) 

Pi = , (5) 



/< = 1 _ ,i- g )g,A^- ((0j )) t/ | (6) 



is the so-called configurational characteristic [|1(J and 



Zt-Y.fr" * ( 7 ) 



stands for the partition function. 

The above procedure, originally employed in overcomes most of the problems posed 
by the old, unnormalized way of evaluating Tsallis'generalized mean values |I|,|TT|j. Some 
hardships remain, though. One of them is that numerical difficulties are sometimes encoun- 
tered, as the Pi expression is explicitly self-referential. An even more serious problem is also 
faced: a maximum is not necessarily guaranteed. Indeed, analyzing the concomitant Hessian 
so as to ascertain just what kind of extreme we face, one encounters the unpleasant fact 
that this Hessian is not diagonal. 

In the present effort we introduce an alternative Lagrange route, that overcomes the 
above mentioned problems. 

II. THE NEW LAGRANGE MULTIPLIERS' SET 

We extremize again (|l]) subject to the constraints (|f), but, and herein lies the central 
idea, rephrase (f|) by recourse to the alternative form 

UI 

i=i 

with j — 1, . . . , n. We have now 

ft \i=i J j=i i=i 
so that, following the customary variational procedure and eliminating chq we find that the 
probabilities are, formally, still given by fl5|). However, in terms of the new set of Lagrange 
multipliers, the configurational characteristics do not depend explicitly on the probabilities 

//=i-(i-ff)EA;(oj°-«o,-»,)- (io) 

3 

Comparing flT0| ) with (^|), it is clear that the Lagrange multipliers Xj of the Tsallis- 
Mendes-Plastino formalism (TMP) and A^- (present treatment) can be connected via 



which leads to the nice result // = fi. The probabilities that appear in ( |TTD are those special 
ones that maximize the entropy, not generic ones. The ensuing, new partition function is 
also of the form (|7]), with fi>0 the well known Tsallis' cut-off condition ||||. Notice that 
now the expression for the MaxEnt probabilities pi is NOT explicitly self-referential. 

In order to ascertain the kind of extreme we are here facing we study the Hessian, that 
now is of diagonal form. The maximum condition simplifies then to the requirement 

d 2 F 

M <0- (12) 
The above derivatives are trivially performed yielding 

H = -«pr 2 /*, (13) 

which formally coincides with the maximum requirement one finds in the case of Tsallis' 
unnormalized formalism. Since the fi are positive-definite quantities, for a maximum one 
should demand that q > 0. 

Extremes found by following the celebrated Lagrange procedure depend only on the 
nature of the constraints, not on the form in which they are expressed. Thus, the two sets 
of multipliers lead to the same numerical values for the micro-state probabilities. Via ( |TTD 
one is always able to establish a connection between both treatments. 

The present algorithm exhibits the same nice properties of the TMP formalism, namely: 

• The MaxEnt probabilities are invariant under uniform shifts of the Hamiltonian's 
energy spectrum (see, for instance, the illuminating discussion of Di Sisto et al. |12|). 

Indeed, after performing the transformation 

€i -»• €i + e (14) 
U q ^U q - e , (15) 

on equation (|5|), with fi given by fllOD, we trivially find that the probabilities pi keep 
their forms invariant if the A^- do not change. Due to relation (|TTD, the Xj are invariant 
too. 



• The mean value of unity equals unity, i.e., = 1, which is not the case with the 
unnormalized expectation values |3|||. 

• One easily finds that, for two independent subsystems A, B, energies add up: U q (A + 
B) = U q (A)+U q (B). 

III. OLM TREATMENT AND RENYI'S MEASURE 

The above OLM recipe has the purpose of diagonalizing the Hessian characterizing the 
kind of extreme one finds after a Lagrange treatment. It is natural to ask whether the 
concomitant recipe works also for other information measures in addition to the Tsallis' 



one. We discuss next the Renyi's measure (R q ) instance |I1]JT3H T5| in an OLM context. We 



have JT6| 



R q _ In [1 + (1 -<?)£,] 
k l-q 

or (Cf. ©) 

R q _ In (Li Pi) 



(16) 



(17) 



k q — 1 

We extremize now R q subject to the constraints with Lagrange multipiers /z*. The 
ensuing Lagrange functional reads 



n w 



p / w \ 



k 

h \i=l / j=l i=l 



whose first derivatives are 

*7 = fe^TES - Q '° - ? "J W - {{ °'^ = a (19) 

The entropy term will lead to a self referencial probability distribution. We have 

(f.*)i^ 

V, = (20) 
with a functional characteristic given by 



/; = ! + (!- ?)£p?£/4 (Of - «Oi» fl ) , (21) 

i 3 

and the partition function 

z; = Y.u:)^- (22) 

i 

We see that the OLM treatment leads to a self-referential solution that the "Renyi-TMP" 
treatament does not exhibit. Indeed, on account of the fact both the constraints and the 
entropy share the J2i Pi denominator, the TMP-Renyi distribution has characteristics of the 
form [15 1 



/i = i + (i-g)E/^(°i <> -« i» s )- ( 23 ) 

j 

Comparing Eqs. (|2"ID and (|23D , one observes that the Lagrange multipliers are connected 



via 



that is the same relation found for Tsallis' measure (see Equation (|TT|)). The following 
correspondence can be appreciated 

Xj = Hj, 

since the configurational characteristics of Equation (p3| ) is identical to the one of Equation 
(0)- 

Consider now, for the second derivatives, the non-diagonal terms. One has. in the OLM 
instance, 

d 2 F* = q 2 { Pl p 3 y- 1 
d Pj d Pi l-q(E kP l) 2 ' 

If one uses the standard TMP approach of Ref. |jj one finds, instead, for these non- 
diagonal terms, 

d 2 F q 2 f l (p lPj )' 1 - 1 



d Pj dpi l-q fckPk) 



9 x2 ' 



(26) 



There is an important distinction to be made between (|25j ) and (|26|), that refers to the 
origin of the non-diagonal terms. In the first case, they have as a source just the entropy 
term of the Lagrangian. In the second, the constraints also contribute. 

As a consequence we gather that the OLM approach will work as nicely as in the Tsallis 
case only for those measures that do not give rise to off-diagonal terms in the associated 
Hessian. 

IV. THERMODYNAMICS 

We pass now to the question of writing down the basic mathematical relationships of 
Thermodynamics, as expressed with respect to the new set of Lagrange multipliers A^. 

In order to do this in the most general quantal fashion we shall work in a basis- 
independent way. This requires consideration of the statistical operator (or density operator) 
p that maximizes Tsallis' entropy, subject to the foreknowledge of M generalized expectation 
values (corresponding to M operators Oj). These take the form 

To these we must add, of course, the normalization requirement 

Trp=l. (28) 

The TMP formalism, where relations are written in terms of the "old" Lagrange multi- 
pliers Xj, yields the usual thermodynamical relationships |IJ, namely 

,, d » (t) = X > ( 29 ) 
!-(ln a Z 5 )=-«A-)),, (30) 

where 



and fllfl 



(32) 



so that the essential mathematical structure of Thermodynamics is preserved. 



Following the standard procedure [0,0] one gets 



P = Zq 



-1 



M 



o, 



1-9 



(33) 



where Z q stands for the partition function 



Z q = Tr 



i-(i-<z)£a;. 6 



o, 



1-9 



(34) 



Enters here Tsallis' cut-off condition [^,TB]. The form ([£J) does not a priori guarantee 
that we will have a positive-definite operator. Some additional considerations are requested. 
Consider the operator 



A =1-(1-?) X>} O 



(35) 



enclosed within parentheses in (|33]). One must ensure its positive-definite character. This 
entails that the eigenvalues of A must be non-negative quantities. This can be achieved by 
recourse to an heuristic cut-off procedure. We replace (|33|) by 



p = z q x \A e(A) 



1/(1-9) 



(36) 



with Z q given by 



Z q = Tr 



A Q(A) 



i/(i-?) 



(37) 



where 9(x) is the Heaviside step-function. Equations (|3Ij|)-(|3"7|) are to be re-interpreted 
as follows. Let |z) and a iy stand, respectively, for the eigenvectors and eigenvalues of the 
operator (B3), whose spectral decomposition is then 



A = ^ai \i) (i\ 



(38) 



In the special basis used above p adopts the appearance 



(39) 



with f(x) defined as 



f(x) = 0, for x<0, 



(40) 



and 



f(x) = x 1 -i , for x > 0. 



(41) 



Notice that f(x) possesses, for < q < 1, a continuous derivative for all x. Moreover, 

df(x) _ ( 1 



\xQ(x)} 



1-9 



dx \q — \ t 

In terms of the statistical operator, Tsallis' entropy S q reads 



(42) 



Q- 


1 


1 




Q- 


1 


z r 


-l 






Q- 


1 



Tr 
Tr 

Tr 



p q {p l - q - i) 
p (z^A e(i) - /) 

Tr(p q ) 



(43) 



p q A Q(A) 



(9-1)' 



where / is the unity operator. 

Obviously, p commutes with A. The product of these two operators can be expressed in 
the common basis that diagonalizes them 

p*i0(i) = V £[/(<*)]* a<|i><*l, (44) 

i 

which entails, passing from the special basis \i) to the general situation, that 



p q AQ{A) = p q 



i-(i-?)E^ ( i-{{° 



(45) 



and, consequently 



Sq -±Tr{p q ) + Z q -^ 3 Tr 

j 



k q 



p q (Oj-UOi 



(46) 



Since the last term of the right-hand-side vanishes, by definition (H), we finally arrive at 

(47) 



q Z^ 1 — 1 
2l - til I Tr (p*) 



k q — 1 

Now, from the very definition (in terms of p) of Tsallis' entropy S q [|5|,|l~6f , we find 



rr(^) = l + (l-?)^, 



so that (g3) and (pD lead to 



and 



(48) 



(49) 



S g = k \n q Z q , 



where \n q Z q has been introduced in (|3lD. 



Using (f£9[), equation ( |TTD can be rewritten as 

A ; 



Following H we define now 



in,z; = in^-^A; O. 



which leads finally to (see (p9|) and (p0|) ) 



<9 



<9 



| ) = = a, 



<9A 



7 (in^;)=-((o,)) c 



(50) 



(51) 



(52) 



(53) 
(54) 



Equations (|5^) and ( |54"D constitute the basic Information Theory relations on which to 
build up, d la Jaynes |7J], Statistical Mechanics. Notar que OLM conduce a las mismas 
relaciones que TMP, y que no ha sido recesario tomar k como constante en ninguna parte 
del desarrollo. 
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As a special instance of Eqs. (p3|) and ([54]) let us discuss the Canonical Ensemble, where 
they adopt the appearance 

d fS„ 



dU q V k 



Z 1 -^' = (55) 



±faz£)=-U q , (56) 



where (see equation (|5^)) 



\n q Z' q = \n q Z q -f3'U q . (57) 



Proponemos que f3' Finally, the specific heat reads 

r _du 1 _ dU^d^ v-adU, ( , R , 

q dT d[3> dT d(3' k'T 2 P d(3 r 1 ' 

We conclude that the mathematical form of the thermodynamic relations is indeed pre- 
served by the present treatment. Both sets of Lagrange multipliers accomplish this feat and 



they are connected via (|TT|). The primed one, however, allows for a simpler treatment, as 
will be illustrated below. 



V. SIMPLE APPLICATIONS 

We consider now some illustrative examples. They are chosen in such a manner that 
each of them discusses a different type of situation: classical and quantal systems, the latter 
in the case of both finite and infinite number of levels. 



A. The classical harmonic oscillator 

Let us consider the classical harmonic oscillator in the canonical ensemble. We can 
associate with the classical oscillator a continuous energy spectrum e(n) = e n with e > 
and n G TZ + compatible with the cut-off condition. The ensuing MaxEnt probabilities adopt 
the appearance 
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where 

f q {n,t') = l-{l-q)^-^-, (60) 

and 

Z q {t')= f nmaX [f q (n,t')] 1/{1 - q) dn, (61) 

J 

with u q = U q /e and t' = k'T/e. We have introduced also n max as the upper integration limit 
on account of Tsallis' cut-off condition. One appreciates the fact that n max — > oo if q > 1. 
nmax is, of course, the maximum n-value that keeps t')] > for g < 1. 

The normalization condition reads 



_ J n — [p g (n,tQ] 9 n^ 

9 Jo^bsMON™ ' 



or, using 



= ^MT^. (63) 

Due to the form of f q , equation fl63f ) constitutes a well-defined expression. By explicitly 
performing the integrals for 1 < q < 2 (for q>2 the integrals diverge) we obtain 

_ t- 2 /(2-g)(l + (l-gK/tQ( 2 -^) 
91 J " t'(l + (1 - q)u q /t')W-<i) ' 1 J 

After a little algebra, the above equation leads to the simple result 

u q {t') = t'. (65) 

Replacing now u q = U q /e and i! = k'T/e, we obtain U q = k'T, so that the specific heat 
reads 

C q = k'. (66) 
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It is worthwhile to remark that, in the case of this particular example, we formally 
regain the usual expressions typical of the q — 1 case. Due to fact that we possess a degree of 
freedom in the definition of k', we can set k' = ks and thus recover Gibbs' Thermodynamics. 
Performing the pertinent integral and using (pop, the partition function becomes 

Z q (t') = t\2-q) l/ ^~ q \ (67) 

According to equation (|o"T|), t 1 can be written in terms of t and Z q , allowing us to recover 

! 

Z q (t) =t l / q (2-q) l l^ 1 '^\ (68) 

and, consequently, 

u q = t 1/q (2-q) 1/q (69) 

C q = |(2-g) 1 ^(i-5)/ g _ ( 70 ) 

These results are identical to those of 0, but are here derived in a remarkably simpler 
fashion. 



B. The two- level system and the quantum harmonic oscillator 

Let us consider the discrete case of a single particle with an energy spectrum given by 
E n = en, where e > and n = 0, 1, ...,N. If N — 1, we are facing the non degenerate 
two level system, while, if n — > oo, the attendant problem is that of the quantum harmonic 
oscillator. 

The micro-state probabilities are of the form, once again 

Pn = J -^~ (71) 

with 

N 

Z q = T,fn m ' 9) - (72) 

n=0 
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The configurational characteristics take the form 



/n(0 



(l-q)(n-u q )/t' 



(73) 



where again (see ( |V A| )), t' = k'T/e and u q = U q /k' . 
Using (|7TD, the mean energy can be written as 

f q/(i-q) ' 

Z^ n =0 Jn 



Un 



(74) 



which, using the explicit form of f n and rearranging terms, allows one to write down the 
following equation 

1 



TV 

E 

n=0 



n — u n 



n - Ur, 



0. 



(75) 



which implicitly defines u q . Notice that one does not arrive to a closed expression. However, 
in order to numerically solve for u q , we just face ([F^). This equation is easily solved by 
recourse to the so-called "seed" methods (cut-off always taken care of), with quick conver- 
gence (seconds). This is to be compared to the TMP instance 0. In their case, one faces a 
non-linear coupled system of equations in order to accomplish the same task. This coupled 



system can be recovered from (|75|) and (71), writing t' in terms of t 



C. Magnetic Systems 

Consider now a very simple magnetic model, discussed, for instance, in [fT7|| : a quantum 
system of N spin 1/2 non- interacting atoms in the presence of a uniform, external magnetic 
field H = Hk (oriented along the unit vector k). Each atom is endowed with a magnetic 

moment Jl = g /x S , where /i = e/(2mc) is Bohr's magneton and S = (h/2) a , with 

^(*) 

a standing for the Pauli matrices. The concomitant interaction energy reads 

N 

n = -J2t ) -H = - 9 -^HS z , (76) 
i=i n 

where S = J2i=i S the total (collective) spin operator. The simultaneous eigenvectors of 
S and 5*2 constitute a basis of the concomitant 2 Ar -dimensional space. We have \S, M), with 
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S = 5,..., N/2, M = —S, ■ ■ ■ ,S, and S = N/2 - [N/2] = (1/2) if N is even (odd). The 
corresponding multiplicities are Y(S, M) = Y(S) = N\(2S + l)/[(N/2 - S)\(N/2 + S + 1)!] 
||17||. We recast the Hamiltonian in the simple form 



n 



with x' = gnoH/3'/h. Our statistical operator can be written as 



l-(l-q)x'(S 



1/(1-9) 



where 



Z g = Tr 



l-{l-q)x'[S, 



S : . 



1/(1-9) 



S K )) ) >0. 



Due to the cut-off condition, 1 — (1 — q)x' lS z - 
The mean value of the spin ^-component is computed according to ( f27| ) 

\\ a '// t - Trip") • 
so that, replacing (|7^) into ( PUD and rearranging then terms we arrive at 

9/(1-9) 



Tr 



1 + (1 - q)x' S, 



s z 



0. 



More explicitly, one has 



N/2 



Ens) J2 

S=8 M=-S L 



1 + (1 - q)x' [M - 
which is the equation to be solved in order to find ({S 2 



9/(1-9) 



M 



S, 



(77) 



(78) 



(79) 



(80) 



(81) 



(82) 



Notice that, once again, one faces just a single equation, that can be easily tackled. If one 



uses instead the TMP prescription (as discussed in [I8"||) one has to solve a coupled, highly 
non-linear system of equations. Such a system can be recovered from ([82]) if one replaces x' 
by x/Tr(p q ) and adds the condition Tr(p q ) from (|78|). 

As in we consider now two asymptotic situations from the present viewpoint. 

For x' — > we Taylor-expand (^) around x' = and find 
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4)>, = ^, (83) 
that leads to an effective particle number 

N° eff = qN, (84) 
as in fl8fl . Following the same mechanism and using ([S3D, one finds that 



Tr(p q ) =2 N ^- q \ (85) 
Remembering that x' = x/Tr(p g ), it is possible to recover the TMP normalized solution 



S z )) q = ^V^ 1 ), (86) 



and 



N° ef f = qN 2 N ^\ (87) 
For x' — > oo, and for < q < 1, expression ( |S2"| ) leads to an equation identical to that of 



E^) E ( M -«^»J = ' (88) 

whose solution reads = -N/2. 



VI. CONCLUSIONS 

In order to obtain the probability distribution pi that maximizes Tsallis' entropy subject 
to appropriate constraints, Tsallis-Mendes-Plastino extremize |I| 

^ = - go Ei* - 1 - E a,- u =l Pl j - ((Oj)) 



vi=l / i=l 



and obtain 
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Pi = -hz — , 

where 

h ^p) 

and Z q is the partition function. Two rather unpleasant facts are thus faced, namely, 

• "pi explicitly depends upon the probability distribution (self- reference) . 

• The Hessian of F is not diagonal. 

In this work we have devised a transformation from the original set of Lagrange multi- 
pliers {Xj} to a new set {A^} such that 

• Self-reference is avoided. 

• The Hessian of F becomes diagonal. 

As a consequence, all calculations, whether analytical or numerical, become much simpler 
than in 0, as illustrated with reference to several simple examples. The primed multipliers 

\' = ^ 
3 Ei P? 

incorporate the in their definition. Since one solves directly for the primed multipliers, 
such a simple step considerably simplifies the TMP treatment. Finally, we remark on the fact 
that the two sets of multipliers lead to thermodynamical relationships that involve identical 
intensive quantities (??). 
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